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Aims. The development of powerful infrared observational technics enables the study of very extincted objects and young embedded star 
forming regions. This is especially interesting in the context of massive stars which form and spend a non negligible fraction of their life still 
jy^ enshrouded in their parental molecular cloud. Spectrophotometric calibrations are thus necessary to constrain the physical properties of heavily 
extincted objects. 

Methods. Here, we derive UBVJHK magnitudes and bolometric corrections from a grid of atmosphere models for O stars. Bessel passbands are 
\^ used. Bolometric corrections (EC) are derived as a function of Teff and are subsequently used to derive BC- spectral type (5 T) and Absolute 
Magnitudes -S T relations. 

Results. Infrared magnitudes and, for the first time, bolometric corrections are given for the full range of spectral types and luminosity classes. 
Infrared colors are essentially constant. (H - K)o is 0.05 mag bluer than previously proposed. Optical calibrations are also provided and are 
I ' similar to previous work, except for (S - V)o which is found to be at minimum -0.28 for standard O stars, slightly larger (0.04 mag) than 
O commonly accepted. 

Conclusions. We present a consistent set of photometric calibrations of optical and infrared magnitudes and bolometric corrections for Galactic 
^ O stars as a function of Teff and spectral type based on non-LTE atmosphere models including winds and line-blanketing. 



^ Key words, stars: fundamental parameters - stars: atmospheres - stars: early-type 

5^ 1. Introduction clusters was proposed by 'Bonne ll. Bate & Zinnecked ( fiQgsl) . 

This key question of the formation of the most massive stars 
Massive stars are known to play an important role in vari- j^^^ triggered a number of observational studies aimed at ob- 
ous fields of astrophysics, from stellar physics to ISM stud- ^^-^-^^ constraints on the properties of the youngest objects 
les and chemical evolution of galaxies, and to cosmological g Crowther & Conti, 2004; Biketal., 2005). Due to the 
issues such as the reionisation of the Universe. In particular, ^^^^^ evolutionary timescale of massive stars, heavily extincted 
the connection between massive stars and star formation is y^^^g ^j^j. forcing regions have to be probed, which requires 
very tight: as a result of their short lifetimes, massive stars jjje use of infrared spectrophotometry, 
are associated with star forming events, and their feedback ef- 
fects (radiation, winds) have a strong impact on star formation Although in principle using only specti-oscopy allows a 
processes. Moreover, their ionising fluxes are responsible for derivation of spectral types and luminosity classes (LC), pho- 
nebular emission lines such as Ly„ or H„, two lines usually tometry can be useful. This is the case when spectra have to 
used to trace star formation dKennicutd Il998t iRusseil et aU be corrected for nebular emission always present in star form- 
l2005h . However, the details of the formation of individual mas- ing regions, rendering the line strength/shape uncertain. As a 
sive stars is still a matter of debate: a standard accretion pro- result, spectral classification and luminosity classes determina- 
cess faces the problem of the strong radiative pressure gener- tions are difficult. Moreover the luminosity is usually derived 
ated by the luminosity of young massive proto-stars, so that from observed magnitudes, extinction and bolometiic correc- 
the mass growth can be stopped at around 10 Mp. Although tions. Estimates of extinction often rely on intrinsic colors of 
progress has be en recently made dYorke & Sonnhalter, 2002: stars while the knowledge of bolometric corrections requires 
iKrumholz et all [2005 ), another scenario in which massive atmosphere models. Hence, accurate intrinsic photometi-y is 
stars form through mergers of low mass protostai's in dense crucial to get access to luminosities. Althoug h such photom- 
etry is usually available in the optical range (" KurucztllQTQt 
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is not the case in the infrared where ca librations are inc om- 
plete. The widely used intrinsic colors of lKoorneed ( ll983l) are 
only given for 06 to 09.5 dwarfs and the latest O supergiants. 
[johnson ( 1966) covers the same range of spectral types / lumi- 
nosity class. 

In this context, the recent development of reliable at- 
mosphere models for massive stars is certainly welcome. 
Indeed, the inclusion of line-blanketing in such models now 
allows realistic prediction of atmospheric structures and emer- 
gent spectra which are used to get qua ntitative constraints 
on the properties of massive stars (ICrowther et all 2002i 
Hillier et al.', '2003* 'Bouret et al.', '2003^, 'Martins et al.', '2004^ 
Rgpolust. Puis & Herrero, 2004; Martins et al., 2005). The grid 
of models computed by Martins. Schaerer & Hilheij fcoOSi) 
(herafter MSH05) and the associated SEDs are especially in- 
teresting since they can be used to compute optical and, most 
importantly, near infrared photometry for the whole range of O 
stars. Together with effective temperature scales, calibrations of 
magnitudes and bolometric corrections as a function of spectral 
type can thus be produced. 

In this paper, we have used the SEDs of MSH05 to calculate 
UBVJHK photometry. In Sect. |2l we present our method and 
gives the results which are discussed in Sect.|3]and summarised 
in Sect.m 

2. Synthetic photometry 

2.1. Method and results 

Synthetic photometry has been computed from the grid of 
atmo sphere models presented by iMartins. Schaerer & HillieJ 
(|2005») ' . From the emergent SED (flux per unit of star surface, 
Fy() we computed the magnitude in each band 



2.5 X log I F^B^dA + constant 



(1) 



with the fil ter passband, according to 

iBessel. CasteUi & Ple3 ( Il998l) . Photometry is thii s com- 
puted in the UBV passbands defined by iBessell ll 19901) 
(Johnson-Cou sins system) and in the JHK passbands of 
iBessel & Brett (.1988J (Johnson-Glass system) for the near-IR. 

Once obtained, these magnitudes were subsequently used 
to determine bolometric corrections for each band from 



bCa = m; 



hoi 



2.5xlog(— ) 



(2) 



where BC^ is the associated bolometric correction, M^' is the 
bolometric magnitude of the sun taken to be equal to 4.75 (rec- 
ommendation of lAUL I1999 I) and L is the luminosity. BC,\ is 
thus computed for each model so that we have bolometric cor- 
rections for the whole range of effective temperatures of O 
stars. These values are shown in Fig.^ we see that, for each 
band, there is a tight coiTelation between BCa and log TeB. A 
simple linear regression of the form 



BCa ^ Ax log(reff) + B 

' Note that the spectral energy 
these models are available at the 
|http://www.mpe.mpg.de/~martins/SED.html_ 



(3) 
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Fig. 1. Bolometric correction in different bands as a function of 
effective temperature for the models of MSH05 (crosses). The 
dashed lines are the regression curves (see Table^for the cor- 
responding parameters). The red dotted line in the plot BCy- 
Teff plot is the relation given by MSH05. 

Table 1. Parameters of the linear regression curves for the BCa 
- Teff relations (A, B, see Eq.|3} together with dispersion cr. 



Band 


A 


B 


cr 


U 


-6.23 


26.46 


0.05 


B 


-6.75 


27.70 


0.04 


V 


-6.89 


28.07 


0.05 


J 


-7.20 


28.85 


0.08 


H 


-7.24 


28.89 


0.09 


K 


-7.24 


28.80 


0.10 



gives the calibration between BCa and effective temperature. 
Parameters A and B of Eq.|3lare given in Table [^together with 
the dispersion cr. The dispersion is larger in the infrared due to 
the increasing contribution of wind emission which introduces 
a scatter (see Sect. 12. 2> . 

The Tetf - scales of MSH05 for various luminosity classes 
were subsequently used to convert effective temperatures into 
spectral types, leading to the calibration of bolometric correc- 
tions as a function of spectral type. The results of this simple 
conversion using the "observational" Teff - scale of MSH05 are 
shown in Fig.|2for dwarf stars. 

Finally, using these new BCa-S T calibrations, the relations 
between log and ST of MSH05 and Eq. |2] we can estimate 
the absolute magnitude in all bands for each spectral type and 
luminosity classes: we obtain the Ma - ST relations (Fig. |2l 
shows the relation for dwarf stars). 

To avoid any confusion, we want to stress that using Eq.^ 
gives directly the absolute magnitudes for each model, and in 
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principle inspection of the spectrum can also provide the spec- 
tral type. In practice, we then have and the spectral type 
for each model. However, we are here interested in producing 
calibrations BC^-ST and -ST for the three main luminos- 
ity classes (dwarfs, giants and supergiants) which are defined 
by specific relations between \ogg and ST (see MSH05). Our 
models do not fall exactly on these relations, some being in 
between two luminosity classes (by construction of the grid). 
Hence, the approach adopted here is aimed at taking this into 
account and is well suited for our purpose, namely producing 
calibrations for each luminosity class. 

The results of the calibrations are gathered in Tables Hand 
|3l In the former, the "observational" Teff scale of MSH05 is 
used, while in the latter, we make use of their "theoretical" Teff 
- ST relation. Inspection of Tables |2] and [Slreveals that adopt- 
ing the "theoretical" T^ff scale of MSH05 changes very little 
the magnitudes: the differences are not larger than 0.03 mag. 
The conclusion is the same for bolometric corrections of the 
earliest spectral types, while later type stars suffer from larger 
differences (up to 0.22 mag). This is expected since late types 
are the ones for which the two effective temperature scales dif- 
fer the most. Note that the colors are basically unaffected by the 
change of effective temperature scale since they depend very 
little on Teff. Since at present it is not clear which T^ff - scale 
better reflects the true properties of O stars, we chose to provide 
both calibrations. 

2.2. Accuracy of calibrations 

In Fig. n and Table [2 we see that the dispersions around the 
average BC - Teff relation increases with wavelength. This is 



a natural consequence of the stronger sensitivity of the SED to 
wind parameters at longer wavelengths. Indeed, massive stars 
are known to emit a significant excess of radiation in the IR- 
radi o range due to free-fre e emission originating in their wind 
(e.g. i Lamers & CassineUii fl999.) . This is illustrated in Fig. |3] 
which shows the variation of the SED when the mass loss rate 
is increased by a factor of 3 (which is the typical uncertainty 
claimed by detailed analysis) in a model for an early super- 
giant. This reveals that any variation in the wind parameters 
will change the wind density, which in turn will affect both the 
continuum level and the strength of near IR emission lines (see 
Fig. 13, leading to a modification of the IR photometry. In or- 
der to estimate the magnitude of this effect, we have run test 
models for a sample of stars for which the mass loss rate was 
increased by a factor of 3. The results are gathered in Table |3 
For such a change in wind density the typical variation in bolo- 
metric correction is ~ 0.02 mag in U, while it can reach 0.2 
mag in K. Note that this is similar to the dispersion of relation 

m 

Fig.[2also shows the calibration of BCy as a function of Teff 
obtained by MSH05. We see that it is slightly different from 
the one presented here, although it is based on the same mod- 
els. The reason for this small discrepancy is the use of a fixed 
effective wavelength (5500 A) in the computation of MSH05, 
while here we use filter curves providing a better measure of 
the flux in the V band. The corresponding new effective wave- 
lengths are usually shorter than 5500 A, leading to a slightly 
smaller magnitude and consequently to a slightly larger bolo- 
metric correction, as seen in Fig.^ Recomputing photometry 
as in MSH05 with the more accurate effective wavelengths of 



Table 2. Optical and infrared magnitudes, colors and bolometric corrections of O stars derived using the "observational" Tgff - scale of MSH05. 



ST 


Mu 


Mb 


Mv 


Mj 


Mm 


Mk 


( U - B)o 


(B-V)q 


( J - H)o 


(H - K)o 




BCb 


SCv 


BCi 


BCh 




03V 


-7.31 


-6.15 


-5.86 


-5.19 


-5.07 


-4.98 


-1.16 


-0.28 


-0.11 


-0.10 


-2.54 


-3.70 


-3.99 


-4.66 


-4.78 


-4.87 


04V 


-7.01 


-5.85 


-5.57 


-4.91 


-4.79 


-4.69 


-1.15 


-0.28 


-0.11 


-0.10 


-2.42 


-3.57 


-3.85 


-4.52 


-4.63 


-4.73 


05V 


-6.69 


-5.54 


-5.27 


-4.60 


-4.49 


-4.39 


-1.14 


-0.28 


-0.11 


-0.10 


-2.29 


-3.43 


-3.71 


-4.37 


-4.48 


-4.58 


05.5V 


-6.56 


-5.42 


-5.14 


-4.48 


-4.37 


-4.27 


-1.14 


-0.28 


-0.11 


-0.10 


-2.22 


-3.36 


-3.64 


-4.29 


-4.40 


-4.50 


06V 


-6.40 


-5.27 


-4.99 


-4.34 


-4.22 


-4.13 


-1.13 


-0.28 


-0.11 


-0.10 


-2.15 


-3.28 


-3.56 


-4.21 


-4.33 


-4.42 


06.5V 


-6.24 


-5.12 


-4.84 


-4.19 


-4.08 


-3.98 


-1.13 


-0.27 


-0.11 


-0.10 


-2.08 


-3.21 


-3.48 


-4.13 


-4.24 


-4.34 


07V 


-6.09 


-4.97 


-4.70 


-4.05 


-3.94 


-3.84 


-1.12 


-0.27 


-0.11 


-0.10 


-2.01 


-3.13 


-3.40 


-4.05 


-4.16 


-4.26 


07.5V 


-5.94 


-4.83 


-4.56 


-3.91 


-3.80 


-3.70 


-1.11 


-0.27 


-0.11 


-0.10 


-1.93 


-3.05 


-3.32 


-3.96 


-4.07 


-4.17 


08V 


-5.79 


-4.68 


-4.41 


-3.78 


-3.67 


-3.57 


-1.11 


-0.27 


-0.11 


-0.10 


-1.86 


-2.97 


-3.24 


-3.87 


-3.98 


-4.08 


08.5V 


-5.62 


-4.52 


-4.25 


-3.62 


-3.51 


-3.41 


-1.10 


-0.27 


-0.11 


-0.10 


-1.78 


-2.88 


-3.15 


-3.78 


-3.89 


-3.99 


09V 


-5.48 


-4.38 


-4.12 


-3.48 


-3.38 


-3.28 


-1.10 


-0.27 


-0.11 


-0.10 


-1.70 


-2.79 


-3.06 


-3.69 


-3.80 


-3.90 


09.5V 


-5.34 


-4.25 


-3.98 


-3.36 


-3.25 


-3.15 


-1.09 


-0.26 


-0.11 


-0.10 


-1.61 


-2.70 


-2.97 


-3.59 


-3.70 


-3.80 


03III 


-7.63 


-6.47 


-6.18 


-5.51 


-5.40 


-5.30 


-1.16 


-0.28 


-0.11 


-0.10 


-2.52 


-3.68 


-3.97 


-4.64 


-4.75 


-4.85 


04III 


-7.49 


-6.33 


-6.05 


-5.39 


-5.27 


-5.18 


-1.15 


-0.28 


-0.11 


-0.10 


-2.39 


-3.54 


-3.82 


-4.49 


-4.60 


-4.70 


05III 


-7.33 


-6.18 


-5.91 


-5.25 


-5.14 


-5.04 


-1.14 


-0.28 


-0.11 


-0.10 


-2.25 


-3.39 


-3.67 


-4.33 


-4.44 


-4.54 


05.5III 


-7.25 


-6.11 


-5.84 


-5.18 


-5.07 


-4.97 


-1.14 


-0.28 


-0.11 


-0.10 


-2.18 


-3.31 


-3.59 


-4.24 


-4.36 


-4.46 


06III 


-7.17 


-6.04 


-5.77 


-5.12 


-5.00 


-4.91 


-1.13 


-0.27 


-0.11 


-0.10 


-2.10 


-3.23 


-3.51 


-4.16 


-4.27 


-4.37 


06.5III 


-7.07 


-5.95 


-5.68 


-5.03 


-4.92 


-4.82 


-1.12 


-0.27 


-0.11 


-0.10 


-2.03 


-3.15 


-3.42 


-4.07 


-4.18 


-4.28 


07III 


-7.00 


-5.88 


-5.61 


-4.97 


-4.86 


-4.76 


-1.12 


-0.27 


-0.11 


-0.10 


-1.95 


-3.07 


-3.34 


-3.98 


-4.09 


-4.19 


07.5III 


-6.93 


-5.82 


-5.55 


-4.91 


-4.80 


-4.70 


-1.11 


-0.27 


-0.11 


-0.10 


-1.87 


-2.98 


-3.25 


-3.89 


-4.00 


-4.10 


08III 


-6.84 


-5.74 


-5.47 


-4.83 


-4.72 


-4.62 


-1.10 


-0.27 


-0.11 


-0.10 


-1.79 


-2.89 


-3.16 


-3.79 


-3.90 


-4.00 


08.5III 


-6.75 


-5.65 


-5.39 


-4.76 


-4.65 


-4.55 


-1.10 


-0.27 


-0.11 


-0.10 


-1.70 


-2.80 


-3.06 


-3.69 


-3.80 


-3.90 


09III 


-6.66 


-5.58 


-5.31 


-4.68 


-4.58 


-4.48 


-1.09 


-0.26 


-0.11 


-0.10 


-1.61 


-2.70 


-2.96 


-3.59 


-3.70 


-3.80 


09.5III 


-6.60 


-5.52 


-5.26 


-4.64 


-4.53 


-4.43 


-1.08 


-0.26 


-0.11 


-0.10 


-1.52 


-2.60 


-2.86 


-3.48 


-3.59 


-3.69 


031 


-7.85 


-6.70 


-6.42 


-5.75 


-5.64 


-5.54 


-1.15 


-0.28 


-0.11 


-0.10 


-2.38 


-3.53 


-3.81 


-4.47 


-4.59 


-4.69 


04I 


-7.82 


-6.68 


-6.40 


-5.74 


-5.63 


-5.53 


-1.14 


-0.28 


-0.11 


-0.10 


-2.26 


-3.40 


-3.68 


-4.34 


-4.45 


-4.55 


051 


-7.79 


-6.66 


-6.39 


-5.73 


-5.62 


-5.52 


-1.13 


-0.28 


-0.11 


-0.10 


-2.13 


-3.26 


-3.54 


-4.19 


-4.30 


-4.40 


05.5I 


-7.78 


-6.66 


-6.38 


-5.73 


-5.62 


-5.52 


-1.13 


-0.27 


-0.11 


-0.10 


-2.07 


-3.19 


-3.47 


-4.12 


-4.23 


-4.33 


06I 


-7.77 


-6.65 


-6.38 


-5.73 


-5.62 


-5.52 


-1.12 


-0.27 


-0.11 


-0.10 


-2.00 


-3.12 


-3.40 


-4.04 


-4.15 


-4.25 


06.51 


-7.76 


-6.65 


-6.38 


-5.74 


-5.62 


-5.53 


-1.11 


-0.27 


-0.11 


-0.10 


-1.94 


-3.05 


-3.32 


-3.96 


-4.08 


-4.17 


071 


-7.76 


-6.65 


-6.38 


-5.74 


-5.63 


-5.53 


-1.11 


-0.27 


-0.11 


-0.10 


-1.87 


-2.98 


-3.24 


-3.89 


-4.00 


-4.09 


07.51 


-7.75 


-6.65 


-6.38 


-5.75 


-5.64 


-5.54 


-1.10 


-0.27 


-0.11 


-0.10 


-1.80 


-2.90 


-3.17 


-3.80 


-3.91 


-4.01 


08I 


-7.73 


-6.63 


-6.36 


-5.73 


-5.62 


-5.52 


-1.10 


-0.27 


-0.11 


-0.10 


-1.72 


-2.82 


-3.09 


-3.72 


-3.83 


-3.93 


08.5I 


-7.73 


-6.64 


-6.37 


-5.74 


-5.63 


-5.53 


-1.09 


-0.26 


-0.11 


-0.10 


-1.65 


-2.74 


-3.00 


-3.63 


-3.74 


-3.84 


091 


-7.70 


-6.62 


-6.36 


-5.73 


-5.62 


-5.52 


-1.08 


-0.26 


-0.11 


-0.10 


-1.57 


-2.66 


-2.92 


-3.54 


-3.65 


-3.75 


09.51 


-7.68 


-6.61 


-6.34 


-5.72 


-5.62 


-5.52 


-1.08 


-0.26 


-0.11 


-0.10 


-1.49 


-2.57 


-2.83 


-3.45 


-3.56 


-3.66 
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Table 3. Optical and infrared magnitudes, colors and bolometric corrections of O stars derived using the "theoretical" Tefr - scale of MSH05. 



ST 


Mu 


Mb 


Mv 


Mj 


Mh 


Mk 


((/ - S)o 


{B-Vh 


U - H)o 


(H - K)o 






BCy 


BC, 


BCh 




03V 


-7.30 


-6.14 


-5.85 


-5.18 


-5.07 


-4.97 


-1.16 


-0.28 


-0.11 


-0.10 


-2.52 


-3.69 


-3.97 


-4.65 


-4.76 


-4.86 


04V 


-7.00 


-5.84 


-5.56 


-4.89 


-4.78 


-4.68 


-1.16 


-0.28 


-0.11 


-0.10 


-2.45 


-3.61 


-3.89 


-4.56 


-4.67 


-4.77 


05V 


-6.69 


-5.55 


-5.27 


-4.60 


-4.49 


-4.39 


-1.15 


-0.28 


-0.11 


-0.10 
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-0.11 


-0.10 


-2.10 


-3.23 


-3.50 


-4.16 
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-4.04 


-4.14 
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-4.48 
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Table 4. Bolometric corrections of models with M increased 
by a factor of 3 (A2,B2,C2,D2) compared to initial mod- 
els (A1,B1,C1,D1). The parameters of the initial models are: 
reff=32210 K, logg=3.26, M^IO"'*''^ Mq yr"', Veo=1960 
km s-i for model Al, reif=33340 K, logg=4.01, M^IQ-^-^^ 
Mo yr"', Voo=2544 km s"' for model Bl, reff=42560 K, 
logg=3.71, M^IQ-^-^^ Mo yr"\ Voo=2538 km s"' for model 
CI, refF=48530 K, logg=4.01, M^IQ-^-^^ Mo yr"', Voo=2977 
km s"' for model Dl. 





BCu 


BCb 


fiCv 


SCj 


BCn 


BCk 


Al 


-1.644 


-2.738 


-2.984 


-3.544 


-3.636 


-3.712 


A2 


-1.654 


-2.750 


-3.006 


-3.519 


-3.588 


-3.607 


Bl 


-1.685 


-2.799 


-3.085 


-3.773 


-3.893 


-4.003 


B2 


-1.685 


-2.800 


-3.087 


-3.776 


-3.892 


-3.999 


CI 


-2.442 


-3.592 


-3.865 


-4.506 


-4.603 


-4.681 


C2 


-2.460 


-3.609 


-3.881 


-4.441 


-4.500 


-4.496 


Dl 


-2.815 


-3.980 


-4.272 


-4.958 


-5.071 


-5.174 


D2 


-2.806 


-3.980 


-4.268 


-4.941 


-5.039 


-5.115 




5000 10^ 1.5x10" 2x10" 2.5x10" 



X[A] 

Fig. 3. Variation of the SED as a function of mass loss rate: 
the black solid curve is for M=1.32 X 10"'^ Mq yr"' while the 
red dotted line is for a mass loss rate three times larger The 
other parameters of the model are reff=48530 K, logg=4.01 
and Vco=2977 km s"' . 



the present study leads to excellent agreement with the magni- 
tudes presented here (differences smaller than 0.01 mag). 

In conclusion, we see that the accuracy of the magnitudes 
and bolometric corrections is usually better than 0. 1 mag and 
reaches 0.2 mag in K band. In any case, this uncertainty is of the 
order of the difference between two spectral sub types within a 
luminosity class. 



3. Comparison with previous studies 

3.1. Absolute magnitudes 

In ter ms of infrared magnitudes, iBlum. Conti & Daminelil 
fcOOO) derived absolute K magnitudes for ZAMS O stars us- 
ing evolii t ionary models CSchaller et al„ .199Z) . V - K from 
Koorneefl (Il983l). and BCv -ST and Jeff ST relations of 
Vacca. Garmanv & Schulll (|l996"). The comparison between 
their results and our relation for dwarfs is displayed in Fig.|3 
We see that our calibration is slightly brighter. Since we cali- 
brated photometry for "normal" O stars, it is not surprising that 
the ZAMS K magnitudes given by Blum. Conti & Daminelij 
(200C| are fai nter, ZAMS stars being less bright than stan- 
dard O dwarfs ('Hanson et alll 19971 'Niemela et al.'.'2006V Note 
that the difference between the two calibrations increases with 
spectral type, which is also what is seen when comparing 
the position of standard O stars to the ZAMS in a HR di- 
agram (e.g. Fig. 12 in MSH05). One could also argue that 
the use of different T^ff - scales is responsible for the ob- 
served differences. However, this is not the case: had we used 
the same approach as Blum. ConJi_&_Damineli (2000) to de- 
rive K magnitudes of ZAMS stars, but using the recent T^ff 
scale of MSH05, we would have found fainter K magni- 
tudes than they find. Indeed, for a given ZAMS evolutionary 
model, i.e. for a given Teff and luminosity, B Cy is similar if 
we use the IVacca. Garmanv & Schulll (11996") BCy - Teff cal- 
ibration or the present one (similar to MSH05 too): the dif- 
ference in bolometric corrections is less than 0.1 mag. Hence 
we would derive a similar V magnitude (within 0. 1 mag) for 
the model. Using V - K from Koorneef (1983) to get the K 
magnitude would then lead to a similar Mk (again within 0. 1 
mag). But now, using the re cent T^^ - scale of MSH05 in- 
stead of the one by Vacca. Garmanv & Schull ( 1996) (as in 
iBlum. Conti & Daminelil l2000l) to convert Teff into spectral 
type gives earlier values: the effective temperatures of MSH05 
are indeed cooler for a given spectral type. This means that us- 
ing th e effective temperature calibration s of MSH05 instead of 
that of lVacca. Garmanv & Sch ul? ("1996) translates to a shift of 
the Mk -5 T relation towards earlier spectral types. In that case, 
the difference with the calibration we provide in the present 
work is even higher Hence, the use of different T^ff - scales is 
not responsible for the shift seen in Fig.l^which is more likely 
attributed to the different evolutionary status in the stars con- 
sidered (ZAMS O stars versus "normal" O stars). 

Calibrations of the V magnitudes as a function of spec- 
tral type have been discussed in MSH05 (see their Sect. 5.1): 
reasonable agreement with previous calibrations (differences 
smaller than 0.4 mag) was found. Inspection of Tables |2 and 
|3l together with Tables 1-6 of MSH05 reveals that the present 
V magnitudes are systematically smaller by 0.02-0.08 magni- 
tudes which is simply due to the better computation of the V 
band photometry in the present study as discussed in Sect. 12.21 
This difference is however well within the typical uncertainty 
of the calibrations. 
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Fig. 4. Absolute K magnitudes as a function of spectral types 
for dwarfs. Filled tri angles are this work, while open triangles 
are the calibration o f iBlum. Conti & Damineli (200 (i). The dif- 
fere nce is mostly due to the fact that Blu m. Conti & DamineUl 
(I2OOO1) give parameters for ZAMS stars while we provide cali- 
bration for "normal" O dwarfs. 

3.2. Intrinsic colors 
3.2.1 . Optical range 

Fig-Elshows the relation (B - V)o - Teff from two sets of atmo- 
sphere models: the CMFGEN ^ models of MSH05 and Kurucz 
models (Kurucz, 1979). On average and for a given T^ff, it 
seems that the new colors are slightly redder, and that the dis- 
persion is larger too. However, this can be explained simply by 
the effect of gravity. As shown by Abbott & Hummer ( 1985j) 
(and confirmed here), the SED of O stars depends on log,g: 
when gravity decreases, the optical flux distributions get red- 
der. This is clearly shown in Fig.|5]where for a given rgft a clear 
sequence of bluer {B - V)q appears when \ogg is increased ^. 
The range of log g covered by the Kurucz models is 3.50-5.00, 
while the grid of MSH05 has 3.20 < logg < 4.25 (see coding 
of symbol size in Fig. |5|i. Thus on average, redder (B - V)o 
are naturally expected for these new models (as seen in Fig.|5}- 
The models of MSH05 also do not reach very blue colors since 
they are restricted to lower gravities compared to the Kurucz 
models. It is also worth noting that the CMFGEN models in- 
clu de wind s, which are know to affect the SED (see Sect. 12.21 
and lOabler et al .l ll989i) . However, the effects are reduced in 
the optical range. A simple inspection of Table |4]i"eveals that 
increasing the mass loss rate by a factor 3 leads to changes in 
B -V (equal to BCv - BCb for a given model) by no more than 

^ For more details on the models, see lHillier & Mi ller ( 1998). 
^ Note however that for the models at low T^g, deviations to this 
general trend appear. 



0. 1 mag. Consequently, the redder colors observed in Fig.|5]for 
the new models can safely be attributed to lower log g. 

Further inspection of Tables |2 and |5]reveals that for stan- 
dard O stars of any spectral type / luminosity class, {B - V)q 
is never bluer than -0.28. This is because normal O stars have 
well defined log g -ST relations for which log g never exceeds 
~ 4.0 (see e.g. Fig. 2 in MSH05). Again, this does not mean 
that very blue (lower than -0.30) {B - V)o colors can not be 
obtained with atmosphere models, but this requires very large 
gravities (log g > 4.25) which are not typical of normal O stars 
(but probably more of very young massive stars). 

How does this result compare to previous studies ? Tlohnsonl 
(1966) derived (B - V)q = -0.30/ - 0.32 for 05-09.5 stars 
fitting the envelope of the position of O stars in a two colors 
diagram (see his Fig. 1). This method assumes that the color is 
almost independent of spectral type (which is indeed the case) 
but Johnson uses a sample of stars covering only part of the 
range of spectral types (especially missing the earliest O stars) 
whi ch may lead to un certainties in the derived average values. 

lFitzgeraldl(ll97d) determined intrinsic colors of O stars also 
using the method of the blue-most envelope in two colors dia- 
gram, but he adopted (B - V\) = -0.32 from Johnson ( 19^. 
This color was determined from studies of O stars in associa- 
tions in which the color excess was determined from photome- 
try of late r type stars. 

Later, IConti. Garmanv & Massev| (Il986h derived colors of 
early type stars in the LMC and found that O stars had (B - 
V)o = -0.30 for most O stars and (B - V)q = -0.24 for late 
supergiants. They also highlighted that the intrinsic colors in 
the LMC were slightly redder than in the Galaxy (while they 
should be similar or slightly bluer due to metallicity effects, see 
below), pointing to a possible problem with the Galactic red- 
dening estimates, and consequently with the Galactic intrinsic 
colors. 

lFitzpatric'i3lll988ll also studied OB supergiants in the LMC 
and derived (B - V)q ^ -0.27 for 03-6 stars and (B - V)o = 
-0.23 for 07-9.7 supergiants, stressing the varia tion of col- 
ors wi thin the different O stars luminosity classes. iFitzpatric^ 
(119881) also compared Stromgren and Johnson intrinsic colors 
in the Galaxy and LMC and highlighted again possible calibra- 
tion uncertainties for the Galaxy. 

We find (B - V)o > -0.28 for Galactic O stars. T his is about 
0.04 mag redder than the early results of Ijohnsoni ([^66) and 
Fitzgerald ( 1970), but in reasonable agreement with the more 
recent study of Fitzpatrick, (.1988). although this study is based 
on LMC stars. We have not computed models at Z = 0.5Zo 
(typical of the LMC), but we do not expect changes in (B - V)q 
larger than ~ 0.01 mag. Indeed, effective temperatures are ex- 
pected to be larger - for a given spectral type - in low metalhc- 
ity en vironments (see e.g. Mokiemet al., 2004: Mass ev et all 
l2004 . but examination of Tables Hand |3] shows that (B - V)o 
is very little sensitive to T^s within the whole range of O stars: 
it changes by no more than 0.02 mag. Hence, (B - V)o should 
be very similar in the LMC and in the Galaxy which is con- 
firmed by the goo d agreement between our results and those of 
lFitzpatrick* ( 1988). Note however that we find (B - V)o slightly 
bluer than Fitzpatrick 11988) for late supergiants. The origin 
of this difference is not clear at present. "Wind effects" can be 
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Fig. 5. (B - V)o as a function of effective tem perature for the 
present work (red open triangles) and from iKurucd ll 19791) 
(black full triangles). For the present work, (B - V)o have been 
computed from the original grid of models and are not taken 
directly from Tables |2] and |3l The size of the symbols scales 
with logg (bigger symbols corresponding to larger gravities). 



excluded since weaker winds - expected at lower metallicity 
- should lead to bluer SEDs and colors (although the effect is 
small in the optical, see above discussion), which is the oppo- 
site of what we see. 

Such a change of the intrinsic B - V colors of standard O 
stars * (at least -0.28 instead of -0.32) has important conse- 
quences for the distance determination of OB associations and 
young clusters. An increase of 0.04 mag in (B - V)o translates 
to a decrease of E{B - V) by the same amount, which then 
implies a reduction of Ay by ~ 0.124 mag (adopting Ry = 
3.1). This means that the distance modulus is reduced by the 
sam e amo unt. For clusters such as Trl6 in the Carina region 
(e.g. lOeGioia et all 1200 li) . this is equivalent to a reduction of 
the distance by 5.5 %. Inversely, if the distance is known, a 
reduction of E{B - V) by 0.04 mag implies a luminosity lower 
by 0.05 dex for a given star. 

Fig-Elshows the comparison between {U - B)o colors de- 
rived in the presen t study as a func t ion of spectral type together 
with the values of lSchmidt- Kaler' ("1982^. Our values are sys- 
tematically redder by ~ 0.05 mag. Fitzgerald ( 1970) derived 
(U - B)o between -1.19 and -1.10 for dwarfs (smaller val- 
ues for early spectral types) and between -1.13 and -1.07 for 
09-9.7 giants and supergiants, in marginal agreement with our 
values. IConti. Garmanv & MassevI lll986l[) found similar val- 



Fig. 6. (U-B) intrinsic colors as a function o f spectral type 
for the present work (solid symbols) and from ISchmidt-KaleJ 
^982). Triangles (squares, circles) are for dwarfs (giants, su- 
pergiants). See text for discussion. 



ues (-1.12 to -1.07) for LMC stars and give (U - B)o=-1.17 
for Galactic stars, arguing again t hat this valu e may be too 
blue (see discussion above). Finall v. Ipitzpatricd (fl9 8 8b derived 
(U - fl)n ^ -1 .04 (resp. -1.08) for 03-6 (resp. 07-9.7) su- 
pergiants. lBessel. CastelU & Plez (.1998) derived (U - B)o col- 
ors for O dwarfs using the ATLAS9 atmosphere models of 
Kurucz ( 19931 11994 . Their values range between -1.22 for 
the earliest O stars to -1.07 for the latest ones (see their Table 
9). This is bluer t han o ur results and similar to the results 
of ISchmidt-Kaleil (|l982'), but can be explained by the differ- 
ent Teff - scales used: adopting the coo ler T^ff - ST re lation 
of MSH05 compared to the relat ion of lCrowtheJ (Il998h used 
bv lBessel. Castelli & Ple z ('l998') translates into a shift of the 
(U - B)o - ST relation towards earlier spectral types. Such a 
shift (of about one spectral sub type) reduces - although not 
comp letely - the diff erence between our relation and those of 
Bessel. Castelli & Plez (1998) and Schmidt-Kaler (19_8^, as 
seen in Fig. |6l The origin of the remaining difference is at- 
tributed to the use of more realistic models in our study. 

On average, om {U-B)q colors are thus consistent with pre- 
vious studies given the improvement in the model atmospheres. 



3.2.2. Near-IR range 



There are less studies of near infrared colors of O stars 
than there are for opt ical colors. The mai n on es are by 
Johnson (1966), Whitteffc" ' ' ' 



We stress again that (B - V)o bluer than -0.28 can be obtained in 
atmosphere models if large gravities (logg > 4.25) are encountered, 
but such gravities are not typical of normal O stars 



van Bred^ lll980l) and iKoorneed 

( 1983). As for optical colors, Johnson ( 196o) used the method 
of t he two colors diagram to derive intrinsic (J - K)o. 
■Whittet & van Breda 1 (Il98fll) proceeded differently and deter- 
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mined intrinsic near IR colors from a sample of 65 stars with 
low extinction (as deduced fr om optical c olors). T heir work 
was done in the Glass system ( lGlassLll974 . Finally. iKoorneei 
( 1983) used a compilation of 127 standard stars in the Southern 
hemisphere to produce near-IR photometry and colors in the 
Johnson-Glass system (J and K in the Johnson system, H in the 
Glass system). 

Fig.0and|8lshow the infrared (J - K)o and (H - K)q colors 
as a function of spectral type for the present work and previous 
studies by Johnson ( 1966), Whittet & van Breda ( 1980) and 
iKoorneeil ( Il983l) . The first interesting point to mention is that 
our study covers the whole range of spectral types and luminos- 
ity classes which was not the case previously. Second, our IR 
colors are almost constant through the range of O stars: we do 
not find any significant variation of near-IR colors with spectral 
type or luminosity class. 

The values of (7 - K)q we find are similar to th e ones by 
[Johnson' ('1966'). This is also true for the colors of iKoorneefl 
( 1983), except for late spectral types for which our colors are 
bl uer (by < 0.1 mag), especial ly for supergiants. The values 
of lWhittet & van Breda I (Il980t) are a little bluer than ours for 
dwarfs, and redder for supergiants. 

As for (H - K)n, there is a reasonably g ood agree ment with 
the va lues of IWhittet & van Breda I lll980l) . whereas 'KoorneeJ 
lll983h predicts colors redder than ours. How ever, the dif- 
ference is rather small (~ O.OSmag). iBlum. Conti & Daminell 
1I2OOOI) studied the stellar content of the HII region W42 and 
give intrinsic H - K colors for ZAMS stars (see their Table 1 ). 
However, for O stars they adopt the value of lKoornee3 ll 1 9 8 3 ) . 

The little differences we observe are likely partly attributed 
to different photometric systems. However, the main conclu- 
sion is that contrary to previous studies, we do not find any 
difference in near-IR colors between dwarfs and supergiants. 

Overall, our results are an improvement over previous the- 
oretical analysis due to the inclusion of line-blanketing and 
winds in the atmosphere models, as well as the use of a better 
Jeff - scale. They present the advantage of covering the whole 
range of spectral types and luminosity classes of O stars in a 
consistent way. They also provide for the first time calibrations 
of bolometric coiTections as a function of spectral type for the 
near infrared range. 

4. Conclusion 

We have derived calibrations of UBVJHK photometry of 
O stars as a function of effective temperature and spec- 
tral type using the recent grid of atmosphere models of 
iMartins. Schaerer & Hillieil bOOS'). UBVJHK photometry was 
comp uted as in [ B essel. Castelli & Plez ( 1998) , us i ng the sys- 
tem of lBessel&BrettI ill 9881) (near IR) and lBessell lll990l) (op- 
tical). 

We provide the first calibrations of near-IR photometry, in- 
cluding bolometric corrections, covering the whole range of 
spectral types and luminosity classes of O stars. Infrared col- 
ors are almost constant. (H - K)q is found to be -0. 10, slightly 
lower (0.05 mag) than the value of lKoornee3 (1 1 9 8 3*). 

Optical photometry is consistent with recent studies. One 
exception is the minimum value of (B - V)o for standard O 
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stars (i.e. with \ogg< 4.0) which is found to be -0.28, slightly 
larger than previously accepted (-0.32). This is important when 
estimating reddening and distances of OB associations since an 
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error of 0.04 mag in color excess amounts to an error of ~ 0. 1 
mag in distance modulus (or ~ 0.05 dex in luminosity). 

These calibrations will be useful to study young massive 
stars embedded in star forming regions and to better understand 
their formation process. 
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